A new exotic state in an old material: A tale of SmB 6 
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We review current theoretical and experimental efforts to identify a novel class ol intermetallic 4/ 
and 5/ orbital materials in which strong interactions between itinerant and predominately localized 
degrees of freedom gives rise to a bulk insulating state at low temperatures, while the surface 
remains metallic. This effect arises due to inversion of even parity conduction bands and odd parity 
very narrow /-electron bands. The number of band inversions is mainly determined by the crystal 
symmetry of a material and the corresponding degeneracy of the hybridized /-electron bands. For 
odd number of band inversions the metallic surface states are chiral and therefore remain robust 
against disorder and time-reversal invariant perturbations. We discuss a number of unresolved 
theoretical issues specific to topological Kondo insulators and outline experimental challenges in 
probing the chiral surface states in these materials. 

PACS numbers: 72.15.Qm, 73.23.-b, 73.63.Kv, 75.20. Hr 
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I. INTRODUCTION 

In the past ten years researchers have been fascinated 
with a peculiar kind of materials: topological insulators 
[THE]. These materials host spin- momentum locked (i.e. 
chiral) metallic surface states, which allow them to re- 
main robust to time-reversal invariant perturbations [7J- 
113] . In addition, these materials, when brought into 
contact with s-wave superconductors, will support Ma- 
jorana fermions that are their own antiparticles |14LH5j . 
The combination of these properties makes topological 
insulators promising platforms for spintronics and quan- 
tum computing applications. At the same time mate- 
rials, which have been proven to possess topologically 
protected metallic surface states, have significant bulk 
conductivity [T6HT9] . Therefore, in this sense they are 
not ideal topological insulators. 

One promising route for discovery of ideal topological 
insulators is to examine materials with strong electron- 
electron interactions. First, the electron correlations may 
fully suppress the bulk conductivity. Secondly, electronic 
interactions may significantly enhance the spin-orbit cou- 
pling, which is responsible for the inversion of the bands 
with opposite parity. In weakly correlated Bi-based topo- 
logical insulators spin-orbit coupling inverts the s- and 
p-bands. In correlated topological insulators we expect 
bands with higher orbital number, either p- and <i-orbitals 
or s, d- and /-orbitals to invert. For example, topological 



Mott insulating state has been theoretically predicted for 
(i-orbital pyrochlore irridates within the extended Hub- 
bard model on a honeycomb lattice [2"DH2"1] , while inver- 
sion between Os d-bands and Ce /-bands leads to topo- 
logical insulator in filled skutterudites |5S] and general 
two-dimensional Kondo system where topological insu- 
lating state is hidden inside the ferromagnetic metallic 
state [26]. It is also worth mentioning the theoretical re- 
alization of various interaction driven topological phases 
in ultracold atom systems and graphene [2"TH2"D] . 

In this article, we will focus on recent theoretical and 
experimental breakthroughs in the search for the ideal 
topological insulator in higher orbital systems. The spe- 
cial attention has been given to already existing /-orbital 
materials [3D], such as CeNiSn, Ce3f3i4Pt 3 , YbBi 2 and 
SmB 6 . These materials, which are called Kondo insula- 
tors, have all the necessary features needed for realizing 
topological behavior: strong spin-orbit coupling, strong 
electron-electron interactions and orbitals with opposite 
parity, Table I. Strong spin-orbit coupling is inherent 
in /-electron systems and it guarantees the inversion 
of the bands at the high symmetry points in the Bril- 
louin zone (BZ). Predominantly localized character of the 
/-electrons furnishes strong Coulomb repulsion between 
them, while hybridization between the even parity con- 
duction electrons and /-electrons leads to the emergence 
of the hybridization gap. Interestingly enough, in some 
Kondo insulators onset of the hybridization gap, observed 



by Raman spectroscopy, has a clear features of the second 
order phase transition. In any case, the opening of the 
hybridization gap does not guarantee an insulating gap, 
of course. However, in Kondo insulators total number of 
conduction and /-electrons per unit cell is even and, con- 
sequently, it immediately follows that Kondo insulators 
are strongly correlated analogues of band insulators. 

TABLE I: Strength of Hubbard interaction U and spin-orbit 
coupling A depending on a type of orbital state 
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U (eV) 1.5 1 1.7 2.1 

a (eV) 0.1 -r 0.2 0.4 -f0.6 0.7vl H2 

This article is organized as follows. In the next Sec- 
tion we will review the theoretical models, which lead to 
the original prediction 31 that Kondo insulators with 
tetragonal or orthorhombic crystalline symmetries can 
naturally become a host to topologically protected metal- 
lic surface states. Section III will be devoted to the review 
of recent experimental and theoretical efforts towards the 
understanding the physics of cubic topological Kondo in- 
sulators and, specifically, SmB 6 . In Section IV we dis- 
cuss the open questions answer to which will deepen our 
understanding of topological Kondo insulators. We sum- 
marize the current status of the field and present our 
conclusions in Section V. 



II. THEORIES OF TOPOLOGICAL KONDO 
INSULATORS 

We Section we review the recent theories of topological 
Kondo insulators. We consider the case when the /-ion 
is in tetragonal crystalline field environment and review 
the theoretical results obtained for this case first. We 
then proceed with the discussion of the theories for cu- 
bic topological Kondo insulators, which are relevant for 
SmB6, YhBi2 and Ce3Bi4Pt3 materials. 

A. tetragonal topological Kondo insulators 

The minimal model of Kondo insulators must involve 
conduction and strongly correlated /-electrons as well as 
hybridization between them. Before we write down the 
corresponding periodic Anderson model, we need to spec- 
ify the /-electron states first. Since most of the tetrago- 
nal systems which are insulating or semi-metallic contain 
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FIG. 1: Large- N mean-field theory results for the topological 
Kondo insulator 43 . Panel (a): Phase diagram calculated 
for the case of doubly degenerate /-ion multiplet, N = 2, as 
a function of the average electronic occupation of the /-level, 
< nf < 1. In the local moment regime, n/ « 1, the weak 
topological insulator (WTI) with the topological invariant 
v — (0; 111) is realized. In the mixed-valence regime, n/ ~ 1, 
Kondo insulator is a strong topological insulator (STI) with 
the topological invariant given by v — (1; 111). When the hy- 
bridization between the conduction and /-electrons becomes 
even stronger, then material is a trivial or band insulator (BI), 
v — (0;000). Panels (b) and (c): phase diagram for TV = 4 
and N — 6 correspondingly. Note that WTI state disappears 
as the the degeneracy of the /-multiplet grows. 

Ce, wc consider the model for the Ce ion in the total an- 
gular momentum J = 5/2 state. The six-fold degenerate 
multiplet is then split into the three Kramers doublets 
with the eigenvectors written conveniently in terms of the 
eigenvectors of the angular momentum projection oper- 
ator J z as follows [32] : 

|/* = ±l) = |±l/2), 

|ju = ±2) = cosa| ± 5/2) - sina| T 3/2), (1) 

|/i = ±3> = sina|±5/2) + cosa| T 3/2). 

where angle a determines the degree of admixture be- 
tween the corresponding orbitals. These eigenvectors can 
be conveniently expressed in terms of the spin part of the 
electron wave function \a &s follows: 
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where A^ na are the Clebch-Gordon coefficients, B^ M are 
some known constants determined from the crystalline 
electric field (CEF) potential. 



The minimal model Hamiltonian for the tetragonal 
Kondo insulator then takes the following form 



# = E e kcLck<x 
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where c, creates an electron in the conduction band 

Kcr 

in a plane-wave state with momentum k, spin a =7, \. 
and energy ek, while /jL creates an /-electron in a state 

with momentum k and the multiplet component fi (TLJ) 

(f) 
and energy e^ . We note that the bandwidth of the /- 

electrons is much smaller than the one for the conduction 
electrons. The third term pi describes the momentum 
dependent hybridization between the conduction and /- 
electrons, while the last terms accounts for strong local 
correlations between the /-electrons on site i. This last 
terms is important as it leads to the local moment for- 
mation. 
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FIG. 2: Two topologically distinct states can be realized in 
our model of tetragonal topological Kondo insulators. The 
first one has a topological invariant v = (1;000) and it cor- 
responds to a strong topological insulator. The second one 
with topological invariant v = (0; 111). 

It is very intuitive that the momentum dependence of 
the hybridization amplitude determines the anisotropy of 
the hybridization gap, which becomes an insulating gap if 
total number of electrons per unit cell is even. Formally, 
the momentum dependence of Vk^o- can be written in 
terms of the spherical harmonic functions [3"2l 151] : 
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where VJ-; are the matrix elements, which can be ex- 
pressed in terms of the corresponding Slater-Koster ma- 
trix elements 1351 . Note that the values of the momentum 



in Q are defined everywhere in the Brillouin zone and 



^T( k ) = \ E *T(R)e ik ' R 



(5) 



R^O 



is a tight-binding generalization of the spherical Harmon- 
ics that preserves the translational symmetry of the hy- 
bridization, Vk = Vk+Gj where G is reciprocal lattice 
vector [31] [34]. Here, R are the positions of the Z near- 
est neighbor sites around the magnetic ion. 

The low energy properties of the model ^ can be an- 
alyzed by the employing the following conjecture: the 
effect of the local correlations between the /-electrons 
leads to the renormalization of the hybridization ampli- 
tude and the shift of the /-energy level: 
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where the renormalization factor Z^ is determined by 
the /-electron self-energy part E w (k, u): 
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Ak = -Tr^Vk] 



(8) 



Then the low-energy model can be diagonalizcd yielding 
a band structure which consists of two doubly degenerate 
bands separated by the momentum dependent energy gap 
A k given by 

1, 

2 

where we suppressed the spin and orbital indices for 
brevity. Since the total number of electrons per unit cell 
is even, we are guaranteed to have the lowest two bands 
to be fully occupied and we have a band insulator pro- 
vided that Ak does not have nodes anywhere in the BZ, 
apart from the high-symmetry points, where it vanishes 
since V_k;/i<r = — Vk. t ^ a - Analysis of the momentum de- 
pendence of the function Ak shows that for /i = ±2, ±3, 
function Ak vanishes at k z — ±k, i.e. gap has two point 
nodes 32 . Therefore, hybridization and insulating gap 
always open for /i = ±1/2, corresponding to a ground 
state multiplet with J z = ±1/2. In passing, we note 
that for systems with lower symmetry (i.e. orthorhom- 
bic), the ground state Kramers doublet will be given by 
linear combination of states ([I]). In this case, despite the 
fact that J z — ±1/2 is nodeless, one can show that hy- 
bridization gap acquires nodes at some finite values of 
momentum k (see discussion in Ref. [32j for details). 
The parity at a high symmetry point can be deter- 



mined by 6 m = sgn(£k* 
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.) [2 E SB [34]. Four in- 



dependent Z<i topological indices (uq; V1V2V3) 



strong (a = 0) and three weak indices (a — 1, 2, 3) can be 
constructed from S m : (i) The strong topological index is 

8 

the product of all eight S m 's: Isti = (~ l)' y ° = ] I 3m = 

?n— 1 

±1; (ii) by setting kj = (where j = x, y, and z), three 
high-symmetry planes, Pj = {k : kj = 0}, are formed 
that contain four high-symmetry points each. The prod- 
uct of the parities at these four points defines the cor- 
responding weak-topological index, /{y TI = ( — l) Va — 
II S m = ±1, a — 1,2,3 with integers correspond- 

ing to the axes x,y and z. The existence of the three 
weak topological indices in 3D is related to a Z2 topo- 
logical index for 2D systems (a weak 3D TI is similar to 
a stack of 2D Z-i topological insulators). Because there 
are three independent ways to stack 2D layers to form a 
3D system, the number of independent weak topological 
indices is also three. A conventional band insulator has 
all of the four indices Isti = ^wti = ^wti = ^wti = +1 
or equivalently (0;000). An index I = (—1) (y a = 1) 
indicates a Z 2 topological state with the odd number of 
surface Dirac modes. In tetragonal KI inversion index 8 m 
of a particular high symmetry point m is negative if the 

conduction band is below the /-band: £k* < £ k v ■ Thus 

(f) 
if £k* = o < e k » _ at the F point, while the remaining 

high symmetry points remain inert, £k^o > e k V, then 
^STi = — 1) and hence the tetragonal Kondo insulator 
is a strong-topological insulator, robust against disorder. 
Fig. [TJa). Weak-topological insulators and topologically 
trivial insulators can in principle be found for different 
band structures and different values of £ k , =0 . 

One can go beyond the phenomenological description 
described above and resort to a more microscopic ap- 
proach. Specifically, we can consider the limit of Uff — > 
00 and work in the restricted phase space by projecting 
out the doubly occupied /-states. Formally, this is ac- 
complished by introducing the slave boson operators |37l - 
[41"] . Then, the mean-field analysis of the resulting model 
can be made by replacing the slave boson operators by 
a number, which must be determined self-consistently. 
The mean-field theory is controlled by the degeneracy of 
the /-orbital multiplet. The large- N mean-field analy- 
sis for the case of tetragonal Kondo insulator have been 
done by Tran, Takimoto and Kim [42| and also by Dzero 
|43j . The results of the slave-boson mean-field theory 
generally agree with the phenomenological approach for 
N = 2, Fig. [ija). In Ref. [13] the mean-field diagram has 
also been obtained for N > 2, Fig. flTb,c). Interestingly, 
it was found that the WTI state is suppressed when the 



degeneracy of the /-multiplet grows, so for N = 6 only 
STI state is present. This result indicates that STI state 
should be preferred for higher symmetry, i.e. cubic, sys- 
tems. We will return to this observation when we discuss 
the cubic topological Kondo insulators below. 




FIG. 3: Dynamical mean field theory phase diagram for the 
two-dimensional (2D) Kondo insulator (adopted from Ref. 
[44] ) . As the strength of the Hubbard interaction Uff be- 
tween the /-electrons increases, there are three topologically 
distinct insulating states. Band insulator (BI) is realized for 
small (compared to the conduction electron bandwidth) val- 
ues of Uff. Topological insulator (TI) with band inversion at 
the r point sets in for intermediate values of Uff . When the 
interaction strength Uff increases further, this state is suc- 
ceeded by a topologically non-trivial 7T-TI state: in this state 
bands of the opposite parity invert at the X points of the 
2D Brillouin zone. At the boundaries between topologically 
distinct states, the system is semi-metallic. 

The phenomenological and mean-field models dis- 
cussed so far assume that the interaction between /- 
electrons is infinitely large. How will the result change if 
we consider the finite values of Uff? Can one study the 
evolution of the topological state as a function of £///? In 
other words, can Kondo insulator be adiabatically con- 
nected to a non-interaction band insulator? This and 
related questions has been recently addressed by Werner 
and Assaad [44J by using the dynamical mean-field theory 
(DMFT) to analyze the two-dimensional Kondo insula- 
tors. 

The physical properties of the model described by the 
Anderson lattice Hamiltonian |3| can be captured, in 
particular, by the single particle propagators. The self- 
energy corrections in the single particle correlation func- 
tions are governed by the Hubbard repulsion between the 



/-electrons. Generally, the /-electron self-energy part 
depends both on momentum and frequency. Within the 
DMFT, however, one considers a single /-electron site 
version of the lattice model pi) lattice, which allows one 
to compute the single particle and higher order corre- 
lation functions exactly [15J . Consequently, the self en- 
ergy parts, which encode the correlation between the /- 
electrons, are all momentum independent. Then, the self 
energy parts for the conduction and /-electrons corre- 
sponding to the lattice problem can be determined self- 
consistently from the solution of the impurity problem. 

Interestingly Werner and Assaad have found j44j that 
increasing the strength of the Hubbard interaction leads 
to the series of transitions between normal insulator 
(small Ujf) into a strong topological insulator with the in 
gap state crossing the Fermi level at the T point and then 
into another strong topological insulator, in which the in- 
gap states cross the Fermi level at the X point, Fig. [3] 
These results clearly indicate that previously developed 
concept of adiabatic connection between uncorrelated 
band insulators and strongly correlated Kondo insulators 
[IS] does not hold at least for the two-dimensional Kondo 
lattice. It remains to be verified, of course, whether this 
conclusion will holds for 3D Kondo insulator. In conclu- 
sion, we note that true tetragonal Kondo insulators still 
remain to be experimentally discovered, since until now 
Ce-based tetragonal " Kondo insulators" , such as CeNiSn 
[17] and CeRuSng [IS], upon improving sample's quality 
become semi-metals. 



B. cubic topological Kondo insulators 

Kondo insulators which consistently show insulating 
behavior in transport measurements - most notably 
SmB 6 , YbB 12 and Ce3Bi4Pt 3 - are all cubic. In this 
section we will review the recent theories of cubic topo- 
logical Kondo insulators, using samarium hexaboride as 
a specific example [15H51] . We note, however, that the 
model we will be discussing below should also hold for 
Ce3Bi4Pt3 although the some parameters can be differ- 
ent. 

The magnetic valence configuration of the Sm ion cor- 
responds to the state with the total angular momentum 
J = 5/2. The six-fold degenerate Sm multiplet is split 
by cubic crystal fields into the Ty Kramers doublet and 
Tg quartet. Consequently the eigenstates of the cubic 



crystalline field Hamiltonian are given by: 
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If T7 is the ground state multiplet the system is a 
scmimetal, since the hybridization gap will have nodes 
in the BZ [35]. Thus, for the insulating state one nec- 
essarily needs Tg to be the ground state multiplet. This 
appears to be indeed the case for SmB 6 (as well as for 
Ce3Bi4Pt3) as evidenced by inelastic neutron scattering 
experiments and Raman spectroscopy [52H55] . 

SmB 6 has a cubic CsCl-like crystal structure, Fig. [4] 
with the B6 clusters located at the center of the unit 
cell. From band-theory calculations [50] [SB] , the Bq clus- 
ters act as spacers which mediate electron hopping be- 
tween Sm sites, but are otherwise inert. In addition to 
band structure results, x-ray photoemission spectroscopy 
(XPS) of SmBg [SI] indicate that the conduction bands 
that hybridize with the localized 4/-orbitals are 5<i-states 
which form electron pockets around the X points. In 
particular, the physics of the 4/ orbitals is governed by 
valence fluctuations involving electrons of the Tg quartet 
and the conduction e g -hole states, 4/ 5 ^ 4/ 6 ± h. The 
conduction states must be d x 2_ y 2 and d^ z 2_ r 2 orbitals of 
e g symmetry. Since the lowest lying state must be a quar- 
tet, one immediately concludes that the band inversion 
happens odd number of times and cubic Kondo insulator 
must be strong topological insulator for a moderate value 
of the hybridization [ST]. Indeed, at the T point both d 
and /-bands remain four-fold degenerate, so that T point 
remains topologically inert. 

Recently, Takimoto [15] has performed the low-energy 
analysis of the Anderson Kondo lattice model ^ prop- 
erly generalized to take into account the realistic band 
structure of SmBg. Specifically, the conduction and /- 
electron spectrum have been derived using the tight- 
binding approximation to the next-nearest-neighbors ap- 
proximation. It is important to emphasize that next- 
nearest-neighbors hopping is needed in order to obtain 
the minimum of the conduction bands at the X points of 
the BZ. The hybridization part of the Hamiltonian has 
also been derived using the tight-binding approximation 
restricted to the nearest-neighbors only. One can justify 
this approximation by noting that the orbital momentum 
AZ = ±1 is transferred from the conduction d-states to 
the /-states in the process of hybridization. The phe- 



nomenological analysis similar to the one of Ref. J3T] of 
the resulting model shows that SmB 6 is a strong topolog- 
ical insulator. An important conclusion one draws from 
the results of Ref. [25] is that as a consequence of the 
band inversion at the X points, there will be three Dirac 
cones on the surfaces perpendicular to the main sym- 
metry axes. Subsequent first-principles study based on 
the local-density-approximation (LDA) plus Gutzwiller 
method |50j has confirmed the results of the phenomeno- 
logical theory [49] . 

Most recently, Alexandrov, Dzero and Coleman (ADC) 
have formulated a general model for the cubic topologi- 
cal Kondo insulators [5T]. The ADC model for the three 
dimensional Kondo insulator can be derived from an ef- 
fective (Uff —> oo) model for the one-dimensional Kondo 
insulator by applying a series of unitary transformations. 
Specifically one considers the quartet of /- and d-holes 
described by and orbital and spin index, denoted by the 
combination A = (a, a) (a = 1,2, a = ±1). Conse- 
quently, the d- and /-states are then described by the 
eight component spinor 



*, 



dxU) 



(10) 



where d\(j) destroys an d-hole at site j, while X \(j) = 
|4/ 6 )(4/ 5 ,A| is the the Hubbard operator that destroys 
an /-hole at site j. The tight-binding Hamiltonian de- 
scribing the hybridized f-d system is then 



ff = £*t(t)/»AA'(Ri-It,)*A'0') 



(11) 



i -J 



in which the nearest hopping matrix has the structure 



h(R) 



h d (R) V(R) 
W(R) fc/(R), 



(12) 



where the diagonal elements describe hopping within the 
d- and /- quartets while the off-diagonal parts describe 
the hybridization between them and R <E (±x, ±y,±z) 
is the vector linking nearest neighbors. The various ma- 
trix elements simplify for hopping along the z-axis, where 
they become orbitally and spin diagonal: 



h l (z)=t l 



1 



Vl, 



V{z) = iV 







(13) 



where I — d, f and rji is the ratio of orbital hopping 
elements. In the above, the overlap between the T s ' 
orbitals, which extend perpendicular to the z-axis is ne- 
glected, since the hybridization is dominated by the over- 

(2) 

lap of the the Y\ orbitals, which extend out along the 



z-axis. The hopping matrix elements in the x and y 
directions are then obtained by rotations in orbital/spin 
space, so that h(x) = U y h(z)U^ and h(y) — U- x h(z)U_ x 
where U y and U- x denote 90° rotations about the y and 
negative x axes, respectively. By construction, the hop- 



ping terms h d ^ (k) in the Hamiltonian (11) could have 



been alternatively obtained from the tight-binding ap- 
proximation within the nearest-neighbor approximation. 
However it is important to keep in mind, that the correct 
band structure for SmBg is still recovered by the proper 
choice of the ratio between the corresponding hopping 
amplitudes [5T] . 

Furthermore, the Fourier transformed hopping matri- 
ces h(k) = X)R^(R) e_lkR can then be written in the 
compact form 



h l (k)=t 



iftfu^+mM^) (i-»»)03(k) 
(i-7B)&(k) MK + mfaO*), 



(14) 
where I — d, f. Here e are the bare energies of the iso- 
lated d and f-quartets, while </>i(k) = c x + c y + 4c 2 , 
02 (k) = 3(c x + c y ) and 3 (k) = y/3(c x - c y ) (c a = 
cos k a , a = x, y, z). The hybridization is given by 



V(k) = 



3{a x + ia y ) V3(cr x - ia y ) 



6 I V3(cx - ivy) &x + i&y 



4cr z 



(15) 



where we denote a a — a a sin k a . Naturally the hy- 
bridization between the even parity d-states and odd- 
parity /-states is an odd parity function of momentum 

v(k) = -n-k). 



The model (111 has been analyzed using the slave- 



boson mean-field approximation and it was shown that 
cubic Kondo insulator is a strong topological insula- 
tor (STI) in complete agreement with previous works 
[191 150] . What is more, ADC have shown that the STI 
state extends well into the local moment regime, simi- 
larly to the results form the large- N theory for tetrago- 
nal Kondo insulators, Fig. 1(c). Lastly, we would like to 
emphasize that the procedure for obtaining the effective 
model for a cubic Kondo insulator [5T] can be easily gen- 
eralized to various type of conduction orbitals and there- 
fore can be used to analyze the topological properties of 
the band structure for other cubic Kondo insulators such 
as YbBi2 and Ce3Bi4Pt3. 



III. EXPERIMENT: SAMARIUM 
HEXABORIDE 

In this Section we present a brief overview of ex- 
periments on the canonical Kondo insulator SmBg [57] . 
There exists a vast amount of experimental literature on 
this material, which merits a separate review paper. Here 
we will focus on discussing the experimental properties of 
importance for possible realization of topologically pro- 
tected chiral surface states. 
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FIG. 4: Crystal structure of SmBe. Sm ions are at the center 
of the unit cell and are surrounded by octahedrons of boron 
ions located at the corners of the unit cell. 

SniBg is the first known and most experimentally stud- 
ied heavy-fermion semiconductor |57j . At high temper- 
atures SmB 6 is a metal with the magnetic susceptibil- 
ity showing the Curie- Weiss like behavior, ^ ~ 1/T 
signaling the existence of the local magnetic moments 
originating from Sm /-electrons. At low temperatures, 
T -C T co h ~ 50 K, SrnBg is a narrow gap semiconductor 
with weakly temperature dependent magnetic suscepti- 
bility [3D1I33]. What is more, the average electronic occu- 
pation of the Sm / states is non-integer and is varying be- 
tween the Sm 2+ (4/ 6 ) and Sm 3+ (4/ 5 5d) configurations. 
This mixed valence state appears as a result of strong 
hybridization between the 5d and 4/ Sm electrons [56] . 
which also leads to formation of the hybridization gap 
Eh f=a 10 meV and many-body insulating gap at E g ~ 5 
meV. 

Most importantly, the saturation of resistivity at tem- 
peratures below T* ~ 5K was consistently observed in- 
dependently by many groups [46] . Resistivity saturation 
was initially interpreted as an extrinsic effect due to for- 
mation of the impurity band. However, almost ten years 
after initial discovery [57] , Allen, Batlogg and Wachter 



[55] on the basis of Hall coefficient data have argued 
that observed values of conductivity were too small for a 
metallic-like system. Moreover, with subsequent accumu- 
lation of the experimental data it became clear that this 
is an intrinsic electronic effect, although the controversy 
regarding intrinsic or extrinsic origin of the effect re- 
mained. This controversy was mainly resolved with pres- 
sure experiments [68 , 69J , which have shown that for pres- 
sure above p* ~ 45 kbar the SmBg recovers its metallic 
properties. Interestingly, it was observed [68] the states 
contributing to low-T conductivity dominate the trans- 
port up to high enough temperatures when the transport 
gap becomes fully suppressed, which means that these 
states must be intrinsic property of a system. Transport 
measurements under pressure and applied magnetic field 
as large as 18 Tesla have shown the negative magnetore- 
sistance [Ap(H)/p(H = 0)] p<p * ~ — H 2 at small pres- 
sure, while for pressures above p* the magnetoresistance 
becomes positive and grows as [Ap(H)/p(H = 0)] p>p * ~ 
_ff 3 / 2 . The activation gap, on the other hand, is very 
weakly dependent on the magnetic field, which indicates, 
albeit indirectly, that the Tg quartet is the ground state 
multiplet for Sm 3+ ions, since it has the smallest value 
of the g- factor, gr 8 =2/7 [70] . 

Since the seminal paper [58] the puzzle of temper- 
ature independent resistivity below T* became clearly 
recognized and its origin remained mysterious for almost 
thirty years. It nevertheless was widely believed that low- 
temperature conductivity in SmBg originates primarily 
from the bulk [6~4Tf6"T] and the surface provides signifi- 
cantly small, if any, contribution to conductivity. From 
our current perspective this is quite remarkable, since the 
first topological insulator could have been discovered al- 
most 30 years ago well before the discovery of the integer 
quantum Hall effect! Nevertheless it took many years 
until very recently to ask the question whether the satu- 
ration of resistivity below 5 K in SmBg is purely surface 
effect. Motivated by the theoretical work [3T], several 
groups independently [59H61] have addressed the issue of 
bulk vs. surface conductivity. The transport [59 , Hall 
effect [60] , tunneling [61] measurements and the angular- 
resolved photoemission spectroscopy [62 unambiguously 
showed that below T* only surface of SmB 6 is conduct- 
ing. What is more, the surface metallic states remains 
surprisingly robust against variations in the surface qual- 
ity. 

Thus, the results of the most recent experiments are 
manifestly in support of the initial prediction [31) that 



Kondo insulators host topologically protected surface 
states. Indeed, robustness of the surface states, the fact 
that their appearance is correlated with the emergence 
of the hybridization gap indicate that SmBg is a strong 
topological insulator [I5H5T] . However, a combination of 
high resolution ARPES and tunneling data are needed 
to directly confirm the chirality of the metallic surface 
states in SmBg. 

IV. OPEN QUESTIONS 

The theoretical and experimental results reviewed by 
us so far bring up a number of issues which need to be 
understood. For example, the Raman spectroscopy data 
[541,155] seem to indicate that opening of the hybridization 
gap in SmBg happens in a way very similar to the second 
order phase transition, i.e. hybridization gap plays a role 
of the mean-field order parameter. Typically mean-field- 
like transitions assume well defined separation of energy 
scales, just like it happens in conventional superconduc- 
tors, for example when the small ratio of the Debye fre- 
quency to the Fermi energy renders the mean-field BCS 
theory to be extremely reliable. In addition, the mean- 
field-like onset of the hybridization implies that the fluc- 
tuations in Kondo insulators are much weaker compared 
with the metallic heavy-fermion systems, which brings 
up the need for better understanding of the fluctuations 
in heavy-fermion systems. 

Theoretically, there is still an open problem of the 
strong coupling description of the topological Kondo in- 
sulators, similar to the Nozieres Fermi liquid picture of 
the Kondo ground state [71]. The progress towards the 
solution of that problem will significantly deepen our un- 
derstanding of the microscopic structure of the chiral 
states on the surface of the topological Kondo insulators. 

Finally, another set of questions questions concerns Ce- 
based heavy fermion semiconductors which order anti- 
ferromagnetically. This materials have tetragonal crystal 
structure and it remains to be seen whether antiferromag- 
netic interactions promote strong topological insulating 
state. The same applies to the possibility for topological 
states in heavy-fermion superconductors. 

V. CONCLUSIONS AND OUTLOOK 

As the search for an ideal topological insulator con- 
tinues, we have come to realizion that at least one ideal 



topological insulator - samarium hexaboride - has been 
discovered almost 30 years ago. Current theories of topo- 
logical Kondo insulators all show that the existence of the 
chiral surface states in /-orbital semiconductors is one of 
its fundamental signatures. Moreover, these states exist 
for the broad range of system's parameters, such as posi- 
tion of the /-electron chemical potential and the strength 
of the Hubbard interaction between the /-electrons. For 
cubic topological Kondo insulators the only requirement 
is that the hybridization between the conduction and 
the /-electrons must be strongest at the X or M high- 
symmetry points in the Brillouin zone. This guarantees 
the odd number of the band inversions and hence a strong 
topological insulator. 

In recent years we have witnessed a resurgence of the- 
oretical and experimental activity in the field of heavy- 
fermions semiconductors and there remains little doubt 
that SmBg is only the first topological Kondo insulators. 
The new potential candidates are YbB 12 and Ce3Bi4Pt3 - 
materials with physical propeties very similar to the ones 
of SmBg. Topological Kondo insulators would present 
an ideal platform for in-depth transport studies of chi- 
ral surface states. In addition, an interplay between 
strong spin-orbit coupling and electron-electron correla- 
tions may open a way to study the broader range of ef- 
fects related to the nontrivial topological structure of the 
electronic states in these materials. One possible direc- 
tion is the in search for topologically non-trivial states in 
Kondo semimetals. However, the most important chal- 
lenge lies in developing new ways to probe the chirality 
of the surface metallic states whether by spin-polarized 
tunneling spectroscopy, by Kerr effect measurements or 
by radio-frequency and microwave spectroscopy. 
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